Feedback Control of Quantum Transport 
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The current through nanostructures hke quantum dots can be stabihzed by a feedback loop 
that continuously adjusts system parameters as a function of the number of tunnelled particles 
n. At large times, the feedback loop freezes the fluctuations of n which leads to highly accurate, 
continuous single particle transfers. For the simplest case of feedback acting simultaneously on all 
system parameters, we show how to reconstruct the original full counting statistics from the frozen 
distribution. 
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Fluctuations of the electronic current have become 
a major tool for probing quantum coherence, inter- 
actions, and dissipation effects in quantum transport 
through nanoscale structures [l|. Monitoring quantum 
objects during their time evolution usually introduces ex- 
tra noise, but it can also compensate backaction effects 
and be used for recycling information in order to control 
the system dynamics J2|. 

The random tunnelling of electrons in quantum trans- 
port is described by the full counting statistics (FCS) 
of transferred charges. Similar to equilibrium thermody- 
namics where, e.g., the cumulants of the particle num- 
ber distribution in the grand canonical ensemble are pro- 
portional to the volume, FCS cumulants in stationary 
transport linearly increase iii^ time (exceptions are pos- 
sible near phase transitions [3]). All quantum transport 
devices thus have to deal with a stochastic element that 
can become a major obstacle when very regular currents 
are required. 

Here, we show that this situation changes by 'freezing' 
the cumulants in time, if one applies feedback (closed 
loop) control [^ to quantum transport. We propose a 
scheme where a time-dependent signal qn{t) is used to 
continuously adjust system parameters such as tunnel 
rates or energy levels. Here, qn{t) = lot — n is an error 
charge determined from the ideal 'target' current Iq and 
the total charge n that has been collected in (or flown out 
of) a reservoir during the measurement (e.g., by a nearby 
quantum point contact detector) up to time t. The error 
charge determines whether to speed up or slow down the 
transport process ~ a form of feedback that is analogous 
to the centrifugal governor used, e.g., in thermo-mechanic 
machines like the steam engine. 

We describe transport in the usual way by coupling 
the system Hamiltonian T-Ls of a few-state nanostructure 
(e.g., a quantum dot) to left and right reservoirs T-Ll/r via 
a tunnel Hamiltonian T-Lt- The feedback loop is modelled 
by an dependence of the parameters in T-Ls and T-Lt on 
time t and on the number operator A'^/j of the right reser- 
voir (drain) . In the usual Born-Markov approximation in 
lowest order in Hr, the reduced system density operator 
p*-"' {t) conditioned on the number of electrons n tunneled 



from left to right (we assume a high-bias situation with 
unidirectional transport) obeys a Master equation 



p(")(t)=£°(t)p(")(t) + J„_i(i)p("-i)(t). 



(1) 



Here, p'"^ {t) is a vector with d real components represent- 
ing system occupations and coherences, and in contrast 
to the usual n-resolved Master equations [J|, the jump 
[J) and non-jump (£°) super-operators {dx d matrices) 
have a time- and n-dependence [5| which can in princi- 
ple be derived from a microscopic model for, e.g., the 
tunnel matrix elements Vk = Vk{t,Nfi) in "Ht, or from 
a corresponding dependence of the energy levels of Hs- 
In an experiment, one could use the signal of a quantum 
point contact in combination with an electronic circuit 
to modulate, e.g., the gate voltages that determine the 
tunnel rates Ti^/ji{n,t) between the reservoirs and the 
nanostructure. 

In all what follows we will assume the elements of £° {t) 
and Sn-i{t) multiplied by analytic functions 



f[qn{t)), qnit) = Iot-n, /(O) = 1 



(2) 



that describe the modulation of system parameters by 
the feedback loop, requiring non-invasive feedback with 
no modulation for zero error charge qn{t) = 0. We also 
assume Iq = I as the stationary current / without feed- 
back in order to simplify some of the formulae; the case 
Iq ^ I yields analogous results. 

Let us start with the simplest transport model: a tun- 
nel junction with no internal system degrees of freedom 
(d = 1) and -£0(i) = Jnit) = r X /(<z„(i)), where F 
is the rate for tunnelling of electrons from left to right. 
We first consider linear feedback f{x) = 1 + gx, where 
5 is a dimensionless feedback parameter. This form is 
appropriate for weak feedback coupling 5 <C 1. A sim- 
ple calculation then yields the first two cumulants of the 
FCS p^"H*) as the average Ci{t) = {n)t = Tt, and the 
variance 



C2it) = {n^)t - {n) 
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(3) 



This already shows that at any finite feedback strength 
g > 0, there occcurs a drastic change: the cumulants 
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FIG. 1: a) Distribution of number n of tunneled particles 
at times t = 30, 60, 100, 140, 180 for tunnel junction (tun- 
nel rate F = 1, linear feedback strength g) ; b) Single re- 
alisation (trajectory) of n/t for tunnel junction as a func- 
tion of time i (F = 1); c) frozen feedback distributions 
Pfb('Ti) = limt^ao pi{n}t + m,t) for chain with TV = 1 and 
N = 2 quantum dots, symbols: numerical solution of Eq. ([T]), 
lines from Eq. ((8l)and re-scaled Eq. ([5]), see text. 



Cfc, fc > 2 no longer increase linearly in time t but con- 
verge to a constant, e.g. C2(oo) = ^. This means that 
the FCS charge distribution no longer spreads out but 
freezes into a stationary distribution with a fixed shape 
that constantly moves to larger n, with a mean value 
{n)t = lot, cf. Fig. dUi). 

In order to obtain p^"'\t) and all the other cumu- 
lants Cfc, we define the Fourier transform p(Xii) = 



The presence of feedback thus transforms the originally 
Poissonian FCS (Cfe = Tt for <? = 0) into a non-diffusive, 
constantly moving distribution at large times. This qual- 
itative change is underlined by the fact that the finite- 
feedback results Eq. dS]), Eq. ^ are non-perturbative in 
the feedback coupling parameter g > 0. We found a fur- 
ther characteristic feature by solving the n-resolved Mas- 
ter equations Eq. ([T]) via the quantum jump method [7|, 
which in the tunnel junction (d = 1) case amounts to a 
simple stochastic algorithm simulating individual exper- 
imental realisations of electron tunneling histories n(i), 
cf. Fig. ([Dd). On long time scales, feedback suppresses 
large deviations of n{t) (we plot n{t)/t for a clearer pic- 
ture) that without feedback lead to the linear increase 
of the cumulants Ck(t) with time. In contrast, on short 
time scales this distinction is barely visible, which is also 
underlined by the fact that the waiting time distribution 
of electron tunneling w{t) [7| reacts much less sensitive 
to feedback (not shown here). 

As a next step, we elucidate the role of internal sys- 
tem degrees of freedom. The simplest case here are spin- 
polarized electrons with sequential tunneling through a 
chain oi N > 1 single resonant levels. We assume strong 
Coulomb blockade, i.e. only one additional electron on 
the chain at a time, and tunnel rates 



F„(n,t)=F„x/(l + .g„(/of-7i)) 



(7) 



that are simultaneously modulated. Results of the nu- 
merical solution of Eq. ([ij are shown in Fig. ([T]:) for 
identical rates Fq = F. The FCS distributions 



lanis (_/fe, we uenne me rourier transiorm pyx,'') = r'^ 

I]n/°^"'H*)e*'^" and use Eq. ^ to derive the partial dif- Pih{m) = lim / 



ferential equation 
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[(n)t+m]XgJ^(x.t) 



(8) 



^^pix,t) = C{x)f(^Iot-^]pix,t), 



(4) 



where £(x) = r(e'^ — 1) and Iq = T (we set the elemen- 
tary charge — e = 1). From the solution of Eq. (|3]), one 
finds the cumulant generating function (CGF) J-"(x, t) for 
linear feedback, 

J^(x,t)^IottX+-Jn{e'^{l 
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Li2{l-e-'>^)], (5) 



where Li2{z) = J^ ^ ln(l — t). From Eq. dS]), we 
find explicit expressions for the FCS via p'^'^\t) ~ 
/-Tr^s""'^^-^^^'*'' and fo^ the cumulants Cu(t) = 
g, ..fc -Fix, ^)lx=o- In the long time limit, one obtains 



Ckit 



= —Bk- 
9 



k>2. 



(6) 



where the Bh 



dx'^ 



are the fc-th Bernoulli- Seki 



x=0 



numbers (the cumulants thus grow rapidly at large k 



are frozen around their maxima that describe moving 
electron number averages {n)t = lot in the large time 
limit. 

Again, we corroborated our numerical results by an 
analytical solution of the feedback master equation for 
the Fourier transformed density operator as in Eq. ((3|), 
where now p{x, t) is a. N + 1-dimensional vector and the 
Liouvillian matrix C{x) has entries — F on the diagonal, 
F on the lower sub-diagonal and Fe*-^ in the right top cor- 
ner, corresponding to counting electrons after they leave 
the N-th. dot in the chain to the right reservoir. At large 
times t, it turns out that C{x) can be replaced by its 
critical eigenvalue Ao(x) with Ao(0) = 0. This approxi- 
mation is in analogy to the long-time limit of the FCS 
for 5 = (no feedback), but for g > it strongly relies 
on the homogeneous coupling of the feedback function to 
all matrix elements in £(%)• Rescaling the tunnel rate as 
F = 7(7V + 1), wc find the CGF J'(x,i) = lnTrsp(x,t) 
as the r.h.s. of Eq. ^ with /q = 7 and the replacements 
g — ;■ g/{N + 1) and x — > x/{N + 1). The feedback-frozen 
distributions p{b{m), Eq. ([5]), obtained in this way are 
in excellent agreement with the numerical results. We 



mention that we carried out this analysis for the hnear 
feedback case that is appropriate for g <^ 1. In general, 
modelling / as an n-th order polynomial leads to n-th 
order PDE systems which for n > 1 make an analytical 
treatment cumbersome. 

The appearance of the eigenvalue Ao(x) in the above 
analysis suggests that for transport through an arbitrary 
system Hs, the original PCS p{n,t) without feedback 
can be reconstructed from the feedback- frozen pn^im) at 
large times t. This can be verified by writing the CGF 
as J-{Xj t ~^ oo) = ixl^i + h{x), leading to 
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(9) 



This, however, is only valid for feedback functions / mul- 
tiplying the Liouvillian C{x) as a whole, cf. Eq. ^ 
('homogeneous feedback'). From Eq. (0), explicit rela- 
tions between the g = Q current cumulants ((/'^)) = 
a(- \k Ao(x)lx=o 8'iid the stationary feedback cumulants 



Ck = 



follow: for weak feedback 



(linear feedback function), the first cumulant 



(10) 



describes the mean target value I^t of the feedback FCS 
at large times t plus an extra charge that flows into the 
drain reservoir due to the mismatch between the target 
current /g and the no-feedback stationary current ((/^)). 
If the latter is known (e.g., by letting the system run 
without feedback for a certain time), Eq. pUj) can be 
used to determine the value of the feedback parameter g 
in an experiment. 

Second, from Eq. ([9]) one recovers the 5 = Fano factor 
F2 = {{l'^)) / m^)) from the frozen g > cumulants as 



F2 = 2gC2 



(11) 



for linear feedback, plus small corrections like 0((7|) for 
feedback functions / with quadratic (522^^) or higher 
terms in their Taylor expansion around x — 0. Corre- 
sponding relations can be easily derived from Eq. ([9]) for 
higher cumulants, e.g., the skewness for linear feedback 
F3^((/3))/((/i))=6g2c|-f3gC3. 

We tested these results and extended them to non- 
homogeneous feedback by considering the single reso- 
nant level model, where at large bias voltage the Mas- 
ter equation reproduces the exact solution for stationary 
transport quantities without feedback. The jump- and 
no-jump operators in Eq. ^ are two-by-two matrices, 
and we parametrized the bare tunnel rates Tj^m between 
dot and left/right electron reservoirs by the asymmetry 
parameter — 1 < a < 1 [81] via F^^ = Tlj^ (trans- 
port from left to right is assumed). Correspondingly, 
we introduced an asymmetry parameter — 1 < 6 < 1 
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FIG. 2: Single quantum dot: (a) second cumulant C2 of frozen 
FCS with right/left tunnel rates Fn = r^i^ and feedback 
couplings gR = Qljt^, 9 = Ql = 0.02; (b) relative change 
of single dot occupation pi; double quantum dot (c) with 
internal coupling Tc = 1, tunnel rates Fl ~ 10, F_r — 1, 
feedback g — 0.01: Fano factor and skewness for homogeneous 
{F2, F3) and inhomogeneous feedback {F2,Fs). 



as gR — gLjTj;, where ga is the dimensionless feed- 
back coupling in the feedback modulated tunnel rates 
r„(7i, t) = F„ X /(I + g^ilot - n), a = L/R. 

In Fig. ([2^) , we tested Eq. (ITT|) as a function of a: 
For homogeneous coupling b — 0, the numerical large- 
time feedback results precisely match the g — Fano 



factor F2 = 5(1- 



For inhomogeneus feedback 6 7^ 0, 



this is no longer the case, as could be expected: the feed- 
back and non-feedback parts of the super-operators in 
Eq. ([1]) are not proportial to each other, and the corre- 
sponding operators in the Fourier transformed equation 
can no longer be replaced by a single eigenvalue Ao(x)- 
The asymmetry parameter b can now be used to fine-tune 
the frozen FCS: for example, negative (positive) b yields 
an overall stronger (weaker) feedback and thus leads to 
sharper (broader) feedback distributions. We found a 
similar scenario for transport through a double quantum 
dot, where we reconstructed the known results for the 
g — Fano factor F2 [9[ and the skewness F3 from the 
numerically obtained frozen homogeneous feedback cu- 
mulants via Eq. (|9]), cf. Fig. ([2j:). In contrast, inhomo- 
geneous feedback that modulates left, right and internal 
tunnel rates (but not the internal bias e) leads to corre- 
sponding expressions ^2 73 with deviations from i^2/3 that 
provide a measure of feedback inhomogeniety. 

In contrast to the FCS, the internal state of the sys- 



tem (i.e., the density operator J2l 



.in) 



(i)) is barely 



changed by the presence of a feedback loop. This can be 
seen in Fig. ([Jt) where we show the relative change of the 
single-dot occupation pi in the stationary limit. For ho- 
mogeneous feedback 6 = 0, this change can be calculated 



analytically by inserting the form p(x, t) 



Jxlot 



rix) 



into Eq. ((4]) and expanding in the feedback coupling g, 
from which equations for the components of r'(0) foUow 
with the resuh 
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Pi{9 = 0) 
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(1 - ^2) , 



(12) 



thus directly relating a (no-feedback) transport quantity 
(the Fano factor F2) with a (feedback) occupation prob- 
ability. For inhomogeneous feedback, the small modula- 
tion of pi again strongly depends on the sign of h. 

The feedback mechanism discussed here leads to non- 
decaying FCS distributions. A potential application 
could be the highly accurate transfer of single electrons at 
minimal errors. An obvious strong competitor of such a 
feedback-controlled device are single-electron pumps (or 
turnstyles, as we are interested in large voltage bias here) 
that are already used for metrological purposes 'lC|. We 
made this comparison more quantitative by modelling a 
single-electron turnstyle as a single level dot with tunnel 
rates rates (a = L,R), r„(i) = laTJ27Li ^(t—taj), with 
an alternating opening of left and right barrier at times 
ta.j such that the time elapsed after k rounds (left-right- 
left) is k xT with T > the period of the turnstyle. We 
compare the turnstyle with the corresponding feedback- 
controlled single-level dot, under the condition that both 
configurations transfer the same average charge per time 
from left to right. At large times, this equivalence leads to 
^ = tanh ( "V ) 7 where we assumed symmetrical tunnel 

rates 7q = 7 and F = Fq, is the bare rate in Eq. ([7]). At 
large times, the turnstyle with its ever broadening FCS 
thus eventually becomes inferior to the (ideal) feedback- 
controlled dot. We define this transition by the time t* 
where the second FCS cumulants C2 of both devices co- 
incide. This condition yields 



n^^icosh^fl^ 
9 V 2 



(13) 



where we used Eq. (fTTj) and -F2 = 5 for the feed- 
back C^^, and the resuh C|"'"P(i) = J.N, N = 

tanh ( -^ ) /(2 cosh [-j-]) that is obtained by consider- 
ing the change of the Fourier-transformed density oper- 
ator p(x) of the turnstyle during one cycle. Interpreting 
iV as the number of electrons transferred through the 
turnstyle before the first counting error occurs, Eq. p^ 
determines the number of electrons n* = Tt* transferred 
through the feedback device before it becomes superior 
to the turnstyle, which for 7T 3> 1 thus happens after 
0{l/g) turnstyle counting errors. 

Finally, we address the question of finite delay times 
that are unavoidable in realistic feedback loops. In gen- 
eral, delay effects in transport will lead to non-Markovian 
feedback Master equations that generalize the r.h.s. of 
Eq. ([ij to integrals and sums over kernels at earlier times 



t' < t and smaller particle numbers n' < n, plus ad- 
ditional inhomogeneous terms ll|. We estimated time- 
delay effects only, using a delay function A(i') = e^* /"^/r 
with delay time t > in the feedback Master equa- 
tion p(")(i) = -F/Jdi'A(i')[/((7n(i - t'))pi^){t - t') - 
[n ^ n — 1)] of the tunnel junction model. Laplace- 
transforming the equations for the first two moments, we 
obtained the first two cumulants in the long time limit, 
Ci{t — > 00) = Vt (which co-incides with the result for 
zero delay r -^ 0), and C2{t -^ 00) = (l/2g) x (1 - 2Fr). 
This indicates that delays that are short on the time-scale 
of the inverse tunneling rate (Ft ^ 1) do not re- install a 
linear increase of the cumulants in time i, but only weakly 
modify the result Eq. (|6]) for the frozen-cumulants. 

To conclude, our results demonstrate that the informa- 
tion contained in the full counting statistics can be frozen 
in by a feedback loop, which can lead to a strong sup- 
pression of fluctuations on long time scales. We expect 
that these predictions can be tested in quantum trans- 
port experiments in the near future. 
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